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We study a Brownian Carnot cycle introduced by T. Schmiedl and U. Seifert [Europhys. Lett. 
81, 20003 (2008)] from a viewpoint of the linear irreversible thermodynamics. By considering the 
entropy production rate of this cycle, we can determine thermodynamic forces and fluxes of the cycle 
and calculate the Onsager coefficients for general protocols, that is, arbitrary schedules to change 
the potential confining the Brownian particle. We show that these Onsager coefficients contain 
the information of the protocol shape and they satisfy the tight-coupling condition irrespective of 
whatever protocol shape we choose. These properties may give an explanation why the Curzon- 
Ahlborn efficiency often appears in the finite-time heat engines. 
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I. INTRODUCTION 



by using the Onsager relations 



Thermodynamics has been developed from the analysis 
of heat engines. Carnot invented an idealized mathemat- 
ical model of heat engines, now called the Carnot cycle, 
and proved that there exists a fundamental upper bound 
of the efficiency of all heat engines, which is given by the 
Carnot efficiency t]q = 1— T c /Th, where Th and T c are the 
temperatures of the hotter and the colder heat reservoirs, 
respectively. But to attain the upper bound, we need 
to operate the heat engines infinitely slowly (quasistatic 
limit) not to cause irreversibility. In the quasistatic limit, 
the heat engine is of no practical use, because the power, 
defined as work output per unit time, becomes 0. 

Practically we need to operate the heat engines in 
a finite time to obtain a finite power. Curzon and 
Ahlborn Q considered a phenomenological model of a 
finite-time Carnot heat engine under the assumptions 
that heat flow obeys the linear Fourier law and that ir- 
reversibility occurs only due to the heat flow (endorc- 
versible approximation) (see also 0, H|), and they de- 
rived that the efficiency at the maximal power rj max of 
that engine becomes 



VGA- 



(1) 



(Curzon- Ahlborn (CA) efficiency) 



Although their derivation of rjQx may seem model- 
specific, they suggested that ry^A closely approximates 
?7max in a few examples of real heat engines [l| and, above 
all, the simple form Eq. (JI} of t]qa implied its universal- 
ity. In fact, subsequent various theoretical studies indeed 
revealed some sort of universality of the CA efficiency 0- 

m. 

Recently Van den Broeck addressed the generality of 
the CA efficiency from a viewpoint of the linear irre- 
versible thermodynamics [l2j ■ He described heat engines 



Ji — L\\X\ 
Ji = L21X1 



L12X2, 
L22X2, 



(2) 
(3) 
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and proved that ??ca is the upper bound of the efficiency 
at the maximal power ?7 max in the linear response regime 
AT -> 0, where AT = T h — T c . This upper bound can be 
attained when the Onsager coefficients T^'s in Eqs. ([2]) 
and © satisfy the tight coupling condition \q\ = 1, where 
q = L12/ \fLv\L22 is called the coupling strength parame- 
ter. Various theoretical models of the heat engines can be 
understood by using the Onsager relations, ranging from 
steady state Brownian motors fl5l - 420j to a macroscopic 
finite-time Carnot cycle [23j]. 

Recently Schmiedl and Seifert suggested an analyti- 
cally tractable model of a Carnot heat engine from which 
useful work can be extracted through the motion of a 
Brownian particle by changing the potential confining 
the particle [23|. In this paper, we call their model a 
Brownian Carnot cycle. By restricting the potential to 
the harmonic form, they analytically showed that the ef- 
ficiency at the maximal power 77 max agrees with t^ca in 
the limit of AT -> 0. (See Sec. II for details.) The re- 
markable feature of their analysis is that they compared 
all protocols when maximizing the power, where we mean 
by protocol the schedule to change the potential. 

Although they showed that their 77 max attains tjca in 
the linear response regime AT — > 0, but the explicit On- 
sager relations Eqs. © and © have not been obtained 
yet, especially how the information of the protocol is re- 
flected in the Onsager relations is unclear. In this paper, 
we calculate the Onsager coefficients Ty 's of the Brown- 
ian Carnot cycle with the harmonic potential for general 
protocols, that is, arbitrary schedules to change the po- 
tential. 

The organization of this paper is as follows. In Sec. II, 
we introduce the Brownian Carnot cycle and review its 
general properties according to [24]. We demonstrate 
the calculations of its Onsager coefficients in Sec. Ill and 
discuss physical implications of them in Sec. IV. We sum- 
marize this study in Sec. V. 
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FIG. 1: Schematic illustration of a Brownian Carnot cycle. A 
Brownian particle immersed in noisy thermal surroundings at 
temperature T r is operated by changing the spring constant 
of the harmonic potential (dashed lines). The shadow regions 
represent the probability distribution p(x, t). (A) and (C) are 
isothermal steps at T r = Th and at T T — T c , respectively. 
(B) and (D) are instantaneous adiabatic steps which do not 
generate irreversibility. During the isothermal step (A) ((C)), 
p(x, t) evolves according to the Fokker-Planck equation from 
Pa (Pb) to p b (p tt ). 



II. MODEL 

Let us introduce the Brownian Carnot cycle [24| . 
Imagine that a one-dimensional overdamped Brownian 
particle is immersed in thermal surroundings at a tem- 
perature T r (r = h, c) and is trapped in a harmonic poten- 
tial. The surroundings behave like noise on the particle 
and fluctuate its position. Therefore its motion becomes 
probabilistic rather than deterministic. In this situation, 
the time evolution of the probability distribution p(x, t) 
of the position x of the particle at time t can be described 
by the Fokker-Planck equation 



dpjx, t) 
dt 



d_ 

^ dx 



f>(,T,f), (4) 



where V(x,t) is the potential energy, which is given by 
V(x, t) = X(t)x 2 /2, using the time-dependent spring con- 
stant X(t). fj, is the mobility of the particle and we set 
the Boltzmann constant fee as unity. We can construct 
a Brownian Carnot cycle in a statistical sense by using 
Eq. Q as follows (see Fig. Q]): (A). Isothermal step at 
Tr = Th (0 < t < t A ): by changing the spring constant, 
we change the distribution function from initial p a (x) to 
Pb(x) for the duration of t A - (B). Adiabatic step: at time 
t = t A , we switch the reservoir at T^ to a colder one at 
T c instantaneously. Because this instantaneous adiabatic 
step gives the probability distribution no time to relax, 



the probability distribution remains p — Pb(x) during this 
step. (C). Isothermal step at T r = T c (t^ < t < tj^ + tc)'- 
by changing the spring constant, we change the distri- 
bution function from pb{x) to p a {x) for the duration of 
tc- (D). Adiabatic step: at time t = tc, we switch the 
reservoir at T c to the one at Th instantaneously again. 
This adiabatic step also keeps the probability distribu- 
tion as p = p a (x) and then the cycle closes. Since X(t) 
determines the schedule to change the potential, we con- 
sider A(t) as the protocol in this Brownian Carnot cycle. 
During the isothermal steps (A) and (C), the probability 
distribution evolves according to the Fokker-Planck equa- 
tion Eq. ((4]) and we also assume that the initial probabil- 
ity p a (x) is Gaussian with the mean and the variance 
w a as p a (x) = exp(— x 2 / (2w a ))/y / 2irw a - Then it can be 
shown that the probability distribution always remains 
Gaussian with the mean in the case of the harmonic 
potential if we initially prepare it so. 

Here we define the internal energy E(t) and the entropy 
S(t) of the particle at time t as 



E(t) = / dx p(x,t)V(x,t), (5) 

J — oo 

/oo 
dx p(x, t) lnp(x, t), (6) 
-oo 

respectively. The entropy S of the particle does not 
change during the adiabatic steps (B) and (D) because 
the probability distribution remains unchanged. Thus 
the adiabatic steps are surely isentropic, which do not 
generate irreversibility during them at all. 

To calculate the work output during the cycle, we 
consider the time evolution equation of the variance 
w(t) = J_ dx p(x, t)x 2 during the isothermal steps (A) 
and (C) by using Eq. ((H) as 



dvM^ 

dt 
dt 



= -2ii\w {A) + 2fiT h (0<t< t A ), (7) 
= -2n\w {c) + 2(j,T c (t A <t<t A + t c ), (8) 



where we have defined the variance during the isothermal 
steps (A) and (C) as and w^ c \ respectively. Then 
the work output in the isothermal step (A) (0 < t < t A ) 
can be calculated by using Eq. ([7]) as 



W {A) 



Ha roo gy 

Jo * 



dt 



i a 
1 

4/i Jo 



dt 2 

t A ( dw (A) \" 
v dt I 



,(A) 



dt 



-Ti[ln w 



(A)]t A 




-W; 



(A) 



(9) 
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where w[ A) = ± J ' A (^^) 2 /w^ dt is the decrease of 
the work by irreversibility, which vanishes in the qua- 
sistatic limit tA — > oo. Assuming that the probabil- 
ity distribution p is always Gaussian with the variance 
w, we have also defined the internal energy change and 
the entropy change during the isothermal step (A) as 
AS {A) = ±{\nw (A X A and AE {A) = ±[\w( A % A , by us- 
ing Eqs ([5]) and ([6]), respectively. The work output in 
the isothermal step (C) (tA < t < tA + tc) can also be 
calculated likewise as 



W {c) 



tA+tc f°» dV 

dt / dx p(x,t) — 

A J — OO °^ 

tA+tc , d\w^ 

dt 

t A dt 2 



1 



t A +t c (- rf^ (c) ■) 



d icTdt+-T c [\nw]Y A +t - 



tA+tc 
tA 



-W; 



(C) 



■ t c as^-ae(°\ 



(10) 



where AS {C) and AE {C) are defined in the same 

way as in the step (A). We need to determine the dy- 
namics of «;W (t) (i = A, C) to calculate and 
explicitly. Since the cyclic change of A(t) asymptotically 
leads to the Gaussian cyclic change of p(x, t), w^(t) 
(i = A,C) is uniquely determined by A(i). Reversely, 
if we determine w^(t) (i = A, C), then the dynamics of 
A(i) is uniquely determined via Eqs. Q and (JS)). There- 
fore we may regard (t) (i — A, C) as the protocol in 
this Brownian Carnot cycle instead of A(t). 

The work outputs and in the adiabatic 

steps (B) (f = tA) and (D) (t = tA + tc) simply be- 
come the change of the internal energy of the system as 
= -AE^ and = -AE^ because the in- 

stantaneous change of the spring constant during the adi- 
abatic steps does not affect the probability distribution. 
Then the total work output W = J2{La during the 



entire cycle (A)-(D) is summed as 



W = -W (A) - W- C) 

irr irr 



(T h - T C )AS, 



(11) 



where we have used relations J2i=A AE^ — and 
Y^iLa AS^ = due to the periodicity of the system 
and have defined AS as 



AS = AS<- A) = -AS ic) 



(12) 



The heat Qh absorbed into the system from the hotter 
reservoir during the step (A) becomes 



Q h = AE^ + = T h AS - W> A >. (13) 
Here we define the power and the efficiency as 

W 



P = W 



tA+tc' 



(14) 



V 



w_ 



(15) 



where the dot denotes a quantity divided by the one-cycle 
period or a quantity per unit time throughout the paper. 
In the Brownian Carnot cycle, the one-cycle period is 

t A + t C - 

After the above setup for the Brownian Carnot cy- 
cle, Schmiedl and Seifcrt considered the maximization of 
the power as follows. Firstly they calculated the opti- 
mal protocol which maximizes the power Eq. (|14[) un- 
der the condition that the durations tA and tc and the 

boundary values w a and Wb are fixed. This optimization 

(A) 

can be realized by minimizing the functionals in Wj rr 

and Wj„ since AS in Eq. (fTTj) is constant. By solving 
Euler-Lagrange equations for these functionals with given 
boundary values w a and Wb, they obtained the optimal 
protocol as 



^(t) =w a (l + Ww b /w a - l)t/t A f 



(16) 



( c \t) = w b (l + Ww a /w b - l)(t - t A )/t c ) 2 -(l7) 



1 A) (c) 

Then, and W-L under this protocol are expressed 



as 



W (A) = i 

irr 

w (c) = i 

irr 



1 

Ta 

1 



(18) 
(19) 



Secondly they considered the further maximization of the 
power of the optimal protocol by changing the durations 
t A and tc- From dP/dtA = dP/dtc = 0, tA and tc can 
be determined as 



tA = tc = 



4(V^ 



M (T h - T C )AS 



Finally, they obtained the maximal power P n 
efficiency at the maximal power ?7 max as 



Pn 



[i \ \a. y/wb/w a j AT 
16{^/wb~- 

m 



2-77C/2' 



(20) 
and the 

(21) 
(22) 



where 77c = 1 — T c /T\y is the Carnot efficiency. This 
result is remarkable because ?y max is independent of the 
boundary values w a and Wb, although P max depends on 
them. Thus, this implies that rj max in Eq. (|22p remains 
unchanged even when the further maximization of the 
power is performed by changing w a and Wb in Eq. (|21j) . 
Therefore, rjmax in Eq. ([22]) gives the efficiency at the 
maximal power under arbitrary protocols of w(t). 

Moreover, by expanding AT, they obtained ?y m ax = 
AT/(2T) + <3(AT 2 ) which is equal to tj C a in the linear 
order of AT. A remark should be added here. The essen- 
tial point of their derivation is to consider the two-step 



4 



maximization of the power: first by the optimal protocol 
for the fixed durations and second by those durations. 
This strategy to find the maximal power is indeed effec- 
tive for the non-linear response regime of the cycle where 
we have no general theory to describe that regime at 
present. But we can use the linear irreversible thermo- 
dynamics in the linear order in AT. As we will show in 
Sec. Ill, we indeed apply the framework of the linear irre- 
versible thermodynamics to this Brownian Carnot cycle 
and can extract more detailed information on the maxi- 
mization of the power in the linear response regime. 



III. CALCULATIONS OF THE ONSAGER 
COEFFICIENTS 

In this Section, we will calculate the Onsager coeffi- 
cients of the Brownian Carnot cycle. Firstly we intro- 
duce the Onsager relations of a general linear irreversible 
heat engine [12j |. Let us consider that the heat engine 
is working under an external force F and a small tem- 
perature difference AT. The work performed against the 
external force F is defined as W = —Fx, where x is 
the thermodynamically conjugate variable of F. In the 
limit of the small temperature difference AT — > 0, we can 
define a thermodynamic force X\ — F/T c ~ F/T where 
T = (Th+T c )/2 and its conjugate flux as J\ — x. We also 
define the inverse temperature difference 1/T C — 1/Tj, ~ 
AT/T 2 as another thermodynamic force X2 and the heat 
flux from the hotter heat reservoir Qh as its conjugate 
flux J2. Moreover, we assume that the Onsager relations 
hold for these fluxes and forces as 133, UBJ] 



J\ — LuXi + L12X2, 
J2 = L21X1 + L22X2, 



(23) 
(24) 



where Lij 's are the Onsager coefficients with the symme- 
try relation L12 = £21- 

To write down the Onsager relations for the Brown- 
ian Carnot cycle, we need to choose the thermodynamic 
fluxes and forces for it. A typical way for choosing them 
is to consider the total entropy production rate a during 
one cycle. Because the system entropy S does not change 
after one cycle as Y^h=a — 0, only the entropy in- 

crease in the reservoirs contributes to a as 



Qh Qc 



Q h w - Qh 



(25) 



In the linear response regime AT — > 0, a can be approx- 
imated as 



-W 



1 



AT 



T (a + l)t A T 2 



Qh, 



(26) 



where we define a = tc/t A a s the ratio of tc to tA and 
we have neglected the higher terms like O(ATW) and 
0(AT 3 Qh), whose reason will be clarified later. The 
expression a = J\X\ + J2X2 by the linear irreversible 



thermodynamics leads to the decomposition of the ther- 
modynamic forces as 



-W _ AT 
, Xo = — — 

rji 1 ^ rp2 



and their conjugate fluxes as 
1 



(a+l)t A 



J2 = Qh- 



(27) 



(28) 



As is clear from these definitions for the thermodynamic 
fluxes and forces, we understand that the higher terms 
like O(ATW) and 0(AT 3 Q h ) do not contribute to the 
entropy production rate a in the linear response regime 
AT -> 0. That's why we neglected them in Eq. ([26]). By 
using these definitions for the thermodynamic fluxes and 
forces, we can write down the Onsager relations of the 
Brownian Carnot cycle as 

1 —W AT 

(a + l)t A = L ^— + L ^YT' ( 29 ) 

-W AT 
Qh =L 21 — + L22j^. (30) 

Now we calculate the Onsager coefficients L\y's. First 
we calculate in and L21 as follows. By using a scaled 
variable s = t/t A (0 < s < I) instead of t (0 < t < 
t A ) during the isothermal step (A), in Eq. (fTTj) 

becomes 



W {A) = -- 



I 1 



1 ( dw w \2 A 

, , / v* 5 ? 2 ' (31) 

t A 4/x 7 w^' t A 



where we call w^ A '(s) = w^ A \t) the protocol shape. In 
Eq. (|3Tj) . we have divided into the part propor- 

tional to the functional of the protocol shape A m and 
the duration t A . W^ TI during the isothermal step (C) in 
Eq. (fTTj) can also be divided as 



irr 



1 I dw (c) \2 n 
y ds I j„ _ °i 



tc 4/i J a w( c ') 



ds=— (32) 
tc y ' 



by using a scaled variable s = (t — t A )/tc (0 < s < 1) 
instead of t (t A < t < t A + t c ) and w; (c) (s) = w^ c \t). 
Then Eq. pT|) can be expressed by using Eqs. (|3"Tj) and 
as 



W 



C'lIT 

a 



ATAS. 



(33) 



By putting AT = in Eq. (|33|). we can obtain the rela- 
tion 



1 



T 



_ —W 

(a + l)t A ~ [Asa + Cirr/a) (a + 1) ' T ' 

which determines the coefficient Lu as 

T 



T11 



(^ irr + C irr /a) (a + 1) ' 



(34) 



(35) 
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Qh at AT = can also be calculated by using Eqs. ([13 
and (03} as 



Qh 



T 2 AS 



-W 



A; 



(A ilT + Cirr/a) (a + 1) T (a + 



r, (36) 



where the second term can be neglected because it is 
0(W 2 ) quantity from Eq. ([34]). Then the coefficient L 2 i 
can be determined as 



T 2 AS 



->i\ 



(Arr + Cirr/a) (a + 1) ' 



(37) 



Next we calculate L12 and L22 likewise. Putting W = 
in Eq. (|33|) . we obtain the relation 



1 



T 2 AS 



AT 

(a + l)t A ~ (A iTI + Cirr/a) (a + 1) " T 5 "' 
which determines the coefficient L12 as 

r 2 AS 



(38) 



(39) 



From Eqs. (|37[) and ([39|) . we find that the Onsager sym- 
metry relation L21 = T12 surely holds as expected. Qh 
at W = thus becomes 



L12 



(A hT + Ci rr /a) (a + 1) ' 



Qh 



T 3 A5 2 



AT 



-4;, 



(Arr + C hT /a) (a + 1) T 2 (a + l)t A s 



(40) 



from Eqs. ([13| and ([38[1 . where the second term can be 
neglected because it is 0(AT 2 ) quantity from Eq. ([38| . 
Then the last coefficient L22 turns out to be 



T 3 AS 2 



(A iTr + Ci rr /a) (a + 1) ' 



(41) 



Note that these Onsager coefficients satisfy the con- 
straints Lu > 0, L22 > and L11L22 ~ L12L21 > 
which come from the positivity of the entropy production 
rate a. We also note that although they are expressed 
in terms of u)'s through A\ n and Ci rr , we could easily 
switch to the A representation by solving the differential 
equation Eqs. and explicitly and substituting the 
solutions w's into A- llr and Ci rr - 



IV. DISCUSSION 

Now that we have calculated the Onsager coefficients 
of the Brownian Carnot cycle in Sec. Ill, we will dis- 
cuss physical implications of them in this Section. To see 
how the Onsager relations are used to describe the effi- 
ciency and the power of the Brownian Carnot cycle, we 
introduce the general theory to describe the heat engines 
governed by the Onsager relations developed in [12| as 
follows. Generally the power P and the efficiency 77 of 
the linear irreversible heat engine are expressed as 



P = W = -Fx = -J1X1T, 
W _ JiX-jT 



n 



(42) 
(43) 



respectively. When X2 determined by the reservoir's 
temperatures and the Onsager coefficients Lij 's are given, 
we can see that only X\ determines the power P and the 
efficiency 77. Since X\ at the maximal power is given by 
dP/dXi = 0, we obtain ry max as 



AT 



?7max — 



where q is defined as 



2T 2 



V Li\L 2 



(44) 



(45) 



which is called the coupling strength parameter. Tak- 
ing into account the restriction — 1 < q < 1 due to 
L\\L 2 2 — L12L21 > which comes from the positivity 
of the entropy production rate a, r7 ma x becomes the up- 
per bound 



Vn 



AT 
2T 



77CA + 0(AT 2 ) 



(46) 



when the tight coupling condition |g| = 1 is satisfied. 

As seen from Eqs. ([3lT) and ([32l , we have shown 
that when we consider a protocol w^(t) by separating 
it into its shape tB™ (s) and its duration U, the On- 
sager coefficients Eqs. (|55|). ([37]). ([59]l and of the 
Brownian Carnot cycle contain the information of the 
protocol shape through A nr and Ci„- and the thermo- 
dynamic flux J\ is the inverse of the one-cycle period 
(a + l)tA = tA + tc- The remarkable feature is that as 
easily confirmed, they satisfy the tight coupling condi- 
tion I q I = 1, independently of the protocol shape, which 
means that the efficiency at the maximal power is always 
the CA efficiency. As we have shown above, this result 
is derived when the power is maximized by changing the 
thermodynamic force X±. However, X± corresponds to 
Ji via Eq. (|29[) at the fixed temperatures and protocol 
shape. Therefore, we can consider that changing X\ is 
equivalent to changing J\ or the one-cycle period. This 
leads to the result that in the present model of the Brow- 
nian Carnot cycle, the efficiency at the maximal power 
is always the CA efficiency when the power is maximized 
by changing only the one-cycle period with the protocol 
shape fixed. If the above feature of the Onsager coef- 
ficients is common to various heat engines, our result 
may explain the appearance of the CA efficiency found 
in [l|, l2lM23j , where the power is maximized by changing 
only the one-cycle period without choosing the optimal 
protocol shape realizing the true maximal power. 

The relation to the study by Schmiedl and Seifert [24| 
reviewed in Sec. II can be understood as follows. Though 
in Eq. (|4l)|). we have seen that the efficiency at the max- 
imal power 77 max under an arbitrary protocol shape is 
t]ca, the maximal power P max should be lower than the 
true maximal power on the space of all protocol shapes. 
Here let us consider the further maximization of P max 
obtained from the Onsager relations. Maximizing the 
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power Eq. (J42j) by changing X\ , the maximal power P n 
is given by 



L 22 AT' 
4T 3 



(47) 



We notice that L22 still depends on cy, -Airr and Ci rl -. 
Then we can further maximize this P max by minimizing 
the functionals in A m and C m with the boundary values 
w a and Wb, and also by changing a, which reduces to an 
inequality 



< 



(ill ^Jw h jw a 



AT 



The equality is realized when a = 1 and A± 



(48) 



Grr — 



/[i, which corresponds to the case of the 



optimal protocol in Eqs. (fT6|) and (fT7|) . The right-hand 
side of Eq. (|4"5)l is the maximal power that the cycle can 
output in the linear response regime at given boundary 
values w a and Wb and is equal to Eq. (|2ip . Thus the result 
in [24| can be reproduced from the Onsager relations. 



its shape and its duration, where we mean by proto- 
col the schedule to change the potential confining the 
Brownian particle. Then, we found that the Onsager co- 
efficients contain the functionals of the protocol shape 
to change the potential and they satisfy the tight cou- 
pling condition, irrespective of whatever protocol shape 
we choose. This result implies that we can attain the 
Curzon-Ahlborn efficiency when maximizing the power 
by changing only the one-cycle period under an arbitrary 
protocol shape. In this sense, the Curzon-Ahlborn effi- 
ciency looks like the Carnot efficiency because the latter 
can also be attained irrespective of whatever protocol 
shape we choose in the quasistatic limit. Although we 
used the harmonic potential to construct the Brownian 
Carnot cycle, we believe that our results could also be 
applied to the models with general potentials and arbi- 
trary dimensions discussed in [24| , where the efficiency at 
the true maximal power agrees with the Curzon-Ahlborn 
efficiency. We expect that our study will stimulate the 
further discussion on the physics of the finite-time heat 
engines. 



SUMMARY 
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